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A model is proposed for assessing grain boundary thermodynamics based on available bulk 
thermodynamic data. A continuous atomic density field and its gradients are introduced to describe the 
grain boundary region. We derive a density-based Gibbs free energy and compute grain boundary phase 
diagram. Grain boundary segregation and the coexistence of the bulk and grain boundary phases are 
discussed. The model is qualitatively demonstrated on the segregation and interfacial phase separation 
in the Pt-Au system. The relationships between the grain boundary’s atomic density, excess free 
volume, and misorientation angle are discussed. The current density-based model is readily applicable 
to study phase equilibria in grain boundaries and opens possibilities to develop micro structure design 
concepts based on segregation engineering. 

Keywords'. Grain Boundary Thermodynamics; Grain Boundary Phase Diagram; Density-based 
Thermodynamics; Grain Boundary Segregation; Segregation Engineering. 
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1. Introduction 

Grain boundaries (GBs) are one of the main sources of heterogeneity in polycrystalline 
microstructures. They interact with various microstructural elements, e.g. GBs’ network [1], secondary- 
phase particles [2, 3], dislocations [4, 5], and vacancies [6, 7] while they themselves possess distinct 
physical properties similar to the bulk phases. GBs can also interact with solute atoms that influence 
both local and global chemistry of the microstructure. Solute segregation to GBs can mediate a whole 
different sort of phenomena such as segregation-assisted GB premelting [8], phase transition [9, 10], 
precipitation [11, 12] and embrittlement [13, 14]. Stabilization of nanocrystalline materials [15, 16] is 
another example where GB segregation can play a critical role. 

The studies on interfacial segregation can be traced back to the works of Gibbs [17], Langmuir 
[18], McLean [19] and Fowler and Guggenheim [20] on the adsorption. Since then, numerous studies 
on the segregation phenomena have been conducted [21-36]. Van der Waals [37], Cahn and Hilliard 
[38] and Cahn [39] worked out free energy functionals that account for heterogeneous interfacial 
features on the mesoscale. Based on this idea, several models for GB segregation were proposed. Ma et 
al [40] developed a model for studying GB segregation in which GB region is described by its reduced 
atomic coordination number. Based on Kobayashi-Warren-Carter’s work [41], Tang et al developed 
models for order-disorder transition [42,43] and phase transition [44] at GBs. Kim et al [45] proposed 
a model for GB segregation using the equal chemical potential condition [46] between the GB and 
adjacent bulk. The common feature of these models is the non-vanishing gradient in GB structure and/or 
composition which is shown to have a great influence on the GB characteristics. Due to this fact, GB 
phases are sometimes referred to as ‘complexions’ [47-49] to emphasize their confined heterogeneous 
nature in contrast to a bulk phase that is -by definition- homogeneous. 


Despite the enormous knowledge accumulated over the last century about the fundamental 
aspects of GBs [35, 36, 50,51], its application at a technical level is awaiting more comprehensive tools 
to allow systematic GB engineering. In particular, the thermodynamic and kinetic databases for GBs 
are lagging behind, when compared to the bulk phases, which is mainly due to the intrinsic complexities 
in the GBs’ nature [50]. Although GB properties are functions of a large space of variables, one may 
not neglect the fact that they are made of the same materials as for the grain interior (bulk), albeit 
deviating in their structure. Hence, parallel to the novel characterization techniques [52, 53] and 
automated simulation methods [54, 55] that are increasingly expanding our ability to study GBs, an 
alternative approach could be to assess the thermodynamic and kinetic properties of the GBs based on 
the thermodynamics and kinetics of their corresponding bulk materials. For this purpose, we need to 
establish a physical framework that allows an approximation of the GB environment with respect to its 
corresponding reference bulk material. In the current study, I propose a continuous atomic density field 
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to describe the GB region. Based on this idea, (i) the GB’s Gibbs free energy is approximated based on 
the bulk thermodynamic properties, (ii) GB segregation and the coexistence of the bulk and GB phases 
are discussed and, (iii) the GB phase diagram is computed. To demonstrate the application of the current 
model, GB segregation and interfacial phase separation in the Pt-Au system are qualitatively studied. 
Furthermore, the relationship between the relative atomic density field and the GB nature in terms of 
GB excess free volume and misorientation angle are discussed. The current density-based model can 
offer a simple practical approach to understand GB phenomena. This model is readily applicable to 
study and design heterogeneous microstructures based on segregation engineering. In a parallel study, 
we exemplified the application of this model in the segregation engineering of Mn in FeMn system 
[56], 

2. A Density-based Model for GB Thermodynamics 

In his seminal work, van der Waals proposed that an interface can be described by a continuous 
density field [37] and its spatial variations. For a GB, that has to accommodate for the geometrical 
mismatch between two adjacent grains, an analogous picture can be portrayed in which the density of 
the GB will be different than that of the bulk materials. This is, on the one hand, even a simpler case 
than the van der Waals’s setup, as the bulk density on the two sides of the GB are the same. On the 
other hand, however, it can be more complex because of the crystalline nature of the bulk material that 
may extend into the GB region. As a first attempt, here we neglect the heterogeneity of GB density (the 
density variation within the GB plane due to the crystallinity) and instead focus on the possibility of 
considering an ‘average’ atomic density parameter and its gradient terms to describe GB region. 

In the following, we derive the density-based Gibbs free energy of a GB in unary (Sec. 2.1) and 
binary (Sec. 2.2) systems. Using the density-based free energy description, we obtain the GB phase 
diagram and segregation isotherm (Sec. 3). The model is demonstrated in the Pt-Au system. Different 
aspects of the current model, in particular the coexistence of the bulk and GB phases and the relationship 
between the GB density, GB excess free volume, and GB misorientation angle, are discussed in Sec. 4. 


2.1 A GB in a pure substance 

The Gibbs free energy functional of a system made of a pure substance A at constant T and p 
can be written as 




— TS a ) df 


( 1 ) 
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where p n (r) is the atomic density field equivalent to the inverse molar volume field H A is 

enthalpy per unit mole with H A = K A + E A + pV A (K A : kinetic energy, E A : potential energy and pV A : 
mechanical energy), and S A is entropy per unit mole, respectively. Van der Waals has shown [37] that 
the potential energy of a heterogeneous system depends on the (local) density as well as (nonlocal) 
density gradients. In this study, we consider a planar GB separating two infinitely large homogeneous 
grains. For this system, breaking the symmetry only normal to the GB plane, the forces can be realized 
considering only one spatial dimension. Thus the potential energy density will be 

1 f x 

E a (x) = E a (- co) + - I f( r ) dr. (2) 

J — 00 

Here E A (—co) is the potential energy inside the homogeneous grain at r = —co and the second term 
describes the energy stored in a given amount of matter when brought from r = — co to a given position 
r — x. f(r) is the sum of all forces acting on point r: 

n 00 

f(r) = I [f(r + q) - f(r - q)]dq (3) 

where ((r ± q) is the force density at point r upon the interaction between the two material points 
separated by a distance ±q. The interaction forces between atoms depend on the interatomic potential 

U[r ± q) with f(r ± q) - dU> '’ g ~ q ^ ■ 

Although the detailed form of the interatomic potentials/forces can vary for different types of 
atoms, it is well-known that the atomistic interactions rapidly decays over the distance between atoms. 
For the sake of our discussion, we consider here a simple functional form 

f(r ± q) = ————— = ip(q) p n [r ± q) (4) 

q 

where a is a material constant, n is the number of atoms, z is a positive exponent and xp(q) - pr . 

Equation (4) is inspired by interatomic force relations -in the absence of the long-range electrostatic 
interactions- with z » 1 [57]. In this formulation, V L is the characteristic coarse-graining volume over 

which the mesoscale atomic density field is measured by p n (r) = -p. Similar coarse-graining 

descriptions are worked out in defining mesoscale free energy formulations [58, 59]. The specific fomi 
of ip(q) can be obtained from atomistic simulations. It is worth noting that second- or higher-order 
dependencies of the interaction forces on the atomic density (including repulsive atomic interactions) 
can also be assumed. In the current treatment, however, we limit ourselves to the simplest first-order 
approximation as proposed in Eq. (4). 
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In the theories of atomistic simulations, the atomic forces are practically calculated up to a cut¬ 
off radius R c beyond which the interatomic forces are neglected [57]. Hence, one reasonable 

4 tc o 

approximation for the coarse-graining volume is V L ~ —Rfi which can relate the atomistic simulation 

length-scale with the current mesoscale density-based formulation. Inserting Eq. (4) in Eq. (3), and 
integrating by part we can write 


n 00 

f(r) = xp(q) [p n (r + q) - p n (r - q)] d q 

Jo 

ip(q) [p n (r + q) - p n (r - q)]]™ - ( xp(q) + ^ ^ 

J o 


r 00 

= - 

J o 


d[p n {r + q) - p n {r - q)] 
dq 


dq 
d q 


d q 


(5) 


with lim i p(q) =0 and lim[p n (r + q) — p n (r — q) ] = 0. Applying the Taylor expansions p n (r ± 

q-> co q-> 0 


q) = p n (r) ± —^q + 

H dr H 2 ! dr 2 


q ± 


3 ! dr 3 


q 3 + ■■■, Eq. (5) gives 


f(r) = 


f 

Jo 


ipiq) 


- ) d _Pn, d ^_Pn 2 
2 dr + dr 3 q + 


_ pO 


r<9p„i 


d 3 Pn 

dr 

-k a 

dr 3 


dq 


(6) 


where we define E® = — 2/ 0 xp(q) d q and k a - f ip[q)q 2 dq . Using Eqs. (2) and (6) with the 
boundary values E A (— co) = E A p n {— co) and (~^rj = 0. we obtain the potential energy as a 

function of atomic density: 


E a (x) = E A p n (x) -y - (7) 

As a consequence of the Tylor expansions, it is clear that the fourth and higher-order spatial 
derivatives of the atomic density field can still contribute to the potential energy. These are, however, 
neglected in the following. The free energy functional in Eq. (1) can be written as 

g A - J' [e° pi + Pn (k a + P v A - ts a ) + y (v Pn ) 2 ] df (8) 
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where we use a three-dimensional notation and / (V p n ) 2 dV — —f p n V 2 p n dV with Vp n = 0 at the 

boundaries of the integral. In Eq. (8) the two terms p n (K A + pV A — TS A ) and E A p 2 x scale differently 

with respect to the atomic density p n . This is an important consequence of the current derivation which 

allows maintaining the continuity of the density profile across the GB. 

The bulk atomic densities on the two sides of the GB can be different, depending on the 

crystallographic planes that meet at the GB plane. This can result in asymmetric GB structure and 

density profile. In the current study, however, we neglect this possibility and assume the same bulk 

density on the two sides of the GB. Thus we can further simplify our description by defining a 

dimensionless density parameter p(x) — pn ^ in which p B = p n (—co). Hence, the general Gibbs free 

Pn 

energy density can be written as 

G a = E b aP 2 + p ( K» + pV A B -TS*)+^-(Vp) 2 (9) 

in which E A = E A (p B ) 2 , k q = K A (p B ) 2 , K A = K A p B , pV A - pV A p B , S A = S A p B and, the Gibbs free 
energy of the homogeneous bulk phase (for p = 1) is 

Ga=E* + k a + pVj? - TS$ = Hi - TS b a . (10) 
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Bulk GB region • Bulk 

(a) 





Figure 1: (a) and (b) Continuous density profile 
across a flat GB with a density p[x = 0) = p GB 
is obtained based on the current model (Eq. (6)). 
The GB width is related to the model parameters 
A a and k a . The current density profile across GB 
is compared versus schematic drawing of order 
parameters from (c) the KWC phase-fieId model 
for GBs [41] usually applied for studying GB 
phase transitions [42-44] and (d) the classical 
phase-field profile (see for instance [60]). (e) The 
energy of a flat GB as a function of its density 
(Eq. (7)) is shown. The solid line indicates the 
likely range of GB densities (see Sec. 4.2 for 
discussions). 


To consider a GB, one can compare a heterogeneous system including a GB against a 
homogeneous bulk system (without a GB). Thus, subtracting (10) from (9) we obtain 

G g a b = E*(p 2 - 1) + (K b +pV B -TS B )(p- 1) + ^ (Vp) 2 . (11) 
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Equation (11) gives an approximation of the GB Gibbs free energy density G gb . For a planar GB 
centered at position x - 0, the equilibrium density profile across the GB will be 

f) 

which satisfies SG ^ = 0 with p(x = 0) = p GB , p(x = ±ij) — 1 and — 0) = 0. Here 2 p = 

2 n J ** is the width of the GB region. Equation (12) gives a continuous atomic density profile across 

the GB. Having a continuous density profile allows deriving the GB free energy such that at least for 
one material point across the GB (at x = 0) the gradient energy term ( Ip) 2 vanishes. Continuity of 

the atomic density profile across the GB is also confirmed by atomistic simulations [56]. Figure 1 
compares the atomic density profile versus previous phase-field models for GBs [41-44] (Fig. 1(c)) as 
well as the classical phase-field profile across an interface [60-62] (Fig. 1 (d)). Further details about Eq. 
(12) are given in Appendix A. 

In the following, p GB , the average relative atomic density within the GB plane, will be referred 
to as ‘GB density’. Using Eqs. (11) and (12) the GB energy y = 2 G% B dr will be 

Y - a 0 (1 - p GB ) 2 (13) 

with a 0 - ^ J—2E a k a . This equation describes the direct relation between y and p GB for an 

equilibrium density profile across the GB. Figure 1 (e) shows the normalized GB energy — as a function 

a 0 

of p GB . The significance of the GB density and its relation with the GB nature is discussed in Sec. 4.2. 


2.2 A GB in a regular binary alloy 

In order to extend the current density-based model for a GB in a binary alloy, we need to discuss 
the significance of mixing energy as a function of atomic density. For a bulk regular solution made of 
A (solvent) and B (solute) atoms, the change in the Gibbs free energy due to the mixing can be 
approximated as [63] 


A GLMb) = A H ? l ix - TAS B nix = H X A X B + RT[X a In + X B In X B ] 

with 

D. - N a Z Ae 


(14) 


(15) 
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in which R is the gas constant, ft is the mixing enthalpy coefficient, N a is Avogadro number and Z is 
the coordination number. Ae = [e A B — 6/1/1 ) 6/iB ) j s bonding energy in which is the bonding energy 
between atoms i and j and X A (X B ) is the mole fraction of atoms A ( B ). In principle, the mixing energy 
terms in Eq. (14) can be assumed to follow the same scaling as found in Eq. (9), i.e. a quadratic p 2 
proportionality for the enthalpy of mixing and a linear proportionality for the entropic contribution. The 
effect of changing the coordination number Z on the GB enthalpy of mixing and GB segregation has 
been previously discussed (see for instance [36, 40, 64] and references therein). While the coordination 
number is proportional to the local density, the stretch in the bonding energy Ae is also expected to 
change with the density in the limit of linear elasticity. Considering these facts and neglecting the 
dependency of configurational entropy for the sake of simplicity, we can write 

A G mix (X BlP ) * p 2 AH B lix - TAS^ ix . (16) 

Following Cahn and Hilliard [38], an energy term due to the concentration gradients -in the presence 
of chemical heterogeneity- enters the free energy description as well. Using Eqs. (9) and (16), the 
density-based Gibbs free energy of a regular solid solution (SS) can be written as 

^ss(^B’P) = XaGa(p) +-^b^b(p) + AG mix (X B ,p) 

- (f> 2 + (Ka+ PVa ~TSa)p + y(Vp) 2 ) 

+ X B (E|p 2 + (K* + pVl - TS B B )p + y (Vp) 2 ) 

+ p 2 nx A x B -TAS^ + ^-m,) 2 

where k x is the concentration gradient coefficient. In Eq. (17), corresponding molar fractions Xj 3 and 
Xf B will be considered for the bulk and GB regions, respectively. For the homogeneous bulk phase, 
one recovers 

g! s (x b b ,p = l) = x a {h a - ts b ) + x b (h b - ts b ) + nx B x B - tas^ lx . (is) 

Using Eq. (17) we are able now to study the GB phase equilibria and segregation isotherms as presented 
in the next section. 


3. Applications of the Current Density-based Model 

3.1 Equilibrium GB phase diagram 
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Obviously, a GB may not exist without its corresponding bulk phase(s). It is, however, useful to 
study the equilibrium phase diagram of a hypothetical GB system as a function of its atomic density. 
As Eq. (17) suggests, for such a hypothetical GB system with p < 1 the phase equilibria should differ 
from that of the corresponding bulk material with p - 1. The equilibrium states are evaluated by 
minimizing the Gibbs free energy functional with respect to the concentration and atomic density fields, 

SQalloy , fiQalloy r* . , SQalloy ^Galloy dQalloy , ~ r ^ _i-» . , 

i.e.-- -» 0 and-- -» 0 with-- - -- — V-- and G nllnv = Gccdr. A numerical 

Sp SX B Sq dq dVq Valloy J SS 

solution of these equations can provide information about the spatiotemporal evolution of phases in the 

GB region. This has been quantitatively demonstrated in a parallel study on the Mn segregation in FeMn 

system [56]. The aim of the current study, however, is to provide analytical descriptions for equilibrium 

GB thermodynamics. Hence, in the following, we rather investigate a single material point at the GB 

center where spatial gradients of the atomic density field vanish, i.e. Vp(x — 0) = 0. We further assume 

that the concentration gradient energy can be neglected. Thus, with p(x — 0) = p GB and X B [x - 0) = 

X B B . the density-based Gibbs free energy description at this material point simplifies to 

g ss {x g b b ,p gb ) = X GB (E B (p GB ) 2 + (K b + pV B - TS B )p GB ) 

+ X% B [E B B {p GB ) 2 + (K B + pv£ - TS B )p GB ) (19) 

+ (p GB ) 2 nx GB x GB -TAS 1 B nix . 

Equation (19) allows to approximate the GB thermodynamic properties, analytically. To 
demonstrate application of the current density-based model, we consider a hypothetical GB with p GB = 
0.75 in the Pt-Au system. In this alloy system, the only difference to the theoretical regular solution is 
in the enthalpy of mixing that extends to the second term in Redlich-Kister polynomial: 

n PtA u = L a + L l (X B t -X B u ). (20) 

Inserting Eq. (20) in Eq. (19) we obtain 

Gss(X™,P GB ) = Xp? ( E B Pt [p GB ) 2 + (K B t + pV B t - TS B t )p GB ) 

+ X G £ (E B u (p GB ) 2 + [K b u + pV B u - TS B u )p GB ) (21) 

+ (p GB ) 2 Cl PtAu X GB X G £-TAS 7 B lix . 
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Figure 2: Gibbs free energy density of the Pt-Au system (Eq. (21)) is shown as a generalized function of 
composition and atomic density at T = 300 K. For p GB = 1 = p B , the Gibbs free energy of the reference bulk is 

recovered. For p GB < 1 the free energy increases. 
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Figure 3: The Gibbs free energy density of the bulk and a hypothetical GB system with p GB = 0.75 in the Pt-Au 
system obtained from Eq. (21) are compared at different temperatures. As expected, the GB Gibbs free energy 

density is higher than the bulk. 


Grolier et al. [65] have assessed the bulk thermodynamic data for binary the Pt-Au and reported 
L 0 = 11625 + 8.3104 T and = —12616 + 5.8186 T [J/mol] for the FCC Pt-Au solid solution. The 
rest of the thermodynamic data for pure Pt and Au are extracted from SGTE compilation by Dinsdale 
[66]. Inserting these values in Eqs. (20) and (21), one obtains the Gibbs free energy as a function of GB 
density and composition. Figure 2 shows a 3D map of the Gibbs free energy for the Pt-Au system at 
300 K as a function of composition Xau and GB density p GB . It is found that as the GB density 
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decreases, the Gibbs free energy increases and deviates larger from the corresponding bulk values. Fig. 
3 compares bulk and GB Gibbs free energy for p GB — 0.75 at different temperatures. 

Using Eq. (21), the equilibrium phase diagrams of the bulk and any specific GB (with a given 
GB density) in the Pt-Au system can be generated. Figure 4 compares the phase diagrams of the bulk 
and a hypothetical GB (with p GB — 0.75) in the Pt-Au solid solution. The results show that, inside a 
(hypothetically standalone) GB, the phase separation becomes possible but for smaller ranges of 
temperature and composition when compared to the bulk phase diagram. 



X Au 


Figure 4: Equilibrium phase diagram for the bulk p = p B = 1 and a hypothetical GB system with p = p GB = 
0.75 are obtained for the Pt-Au system. The dash lines represent the chemical spinodals. The colored area marks 
the possible domain where GB phase diagrams for 0.75 < p GB < 1 can appear. Note that the equilibrium bulk 
and GB phase diagrams are plotted independently. The coexistence of the bulk and GB phases will be discussed 

in Secs. 3.2 and 4.1 and Fig. 6. 


As it is evident from the Gibbs free energy plots (Figs. 2 and 3), the GB thermodynamics strongly 
depend on the GB density p GB . Although the phase diagram is a nonlinear function of the GB density, 
the possible values for GB density p GB are finite and vary close to 1. For instance, the colored area in 
Fig. 4 indicates where the GB phase diagram can appear for 0.75 < p GB < 1. Thus, once the range of 
GB density is determined, the equilibrium GB phase diagrams can be approximated. The GB density 
represents the nature (type) of the GB and correlates with the GB energy and misorientation angle. 
These aspects of the current density-based model are discussed in Sec. 4.2. 
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3.2 GB segregation behavior 

In Section 3.1 a hypothetical standalone GB system was discussed. For a complete description 
of GB phase equilibria, however, the coexistence between the GB and bulk phases must be understood. 
In a binary system, this requires equality of the (relative) chemical potentials pf — Ea = Mb B — Ma B a U 
across the system (parallel tangent construction) [20, 67]. In reality, these are not sufficient conditions 
as a minimum energy state with respect to the atomic density field must be satisfied as well. However, 
if we consider the same material point at the GB center as above, we can simplify the problem and solve 
for the equality of the (relative) chemical potentials using Eq. (19) which gives 


X G b B 

1-X GB 


■xl 


■exp 


[A E b + n] ( p GB2 - l) + (A K b + pAV B - TAS B )(p GB - 1) + 20 [. 


Xl 


-p gb x g b b \ 


RT 


( 22 ) 


Here X% B is the composition of the GB with p = p GB . Xji is the composition of the bulk far from the 
GB, A E b = El - E b a , A K b = K% - K% , pAV B = pE# - pV B , and AS = Sjj - S% . Equation (22) 
resembles the Fowler-Guuggenheim segregation isotherm [20] but also takes the specific effect of a 
given GB, represented by its atomic density, into account. For the Pt-Au system, we obtain 


Xl 


yGB 
' a Au 


■ y B 

A Au 


(23) 


■ exp - 


[A E b + L 0 ] ( p GB2 - l) + (A K b + pAV B - TAS B )(p GB - 1) + 2L 0 [x% u 


?] + <?' 


RT 


with Q = L X [2 P gb 2 X gb (1 - X GB ) - (1 - 1X GB ) 2 - 2X B U ( 1 - X B J + (1 - 2X B U ) 2 ] , A E B = 

e au ~ E Pt » &K B = Kl u - K$ t , pAV B = pV B u - pV$ t , and AS = 5| u - Sf t . Equations (22) and (23) 
use GB density and bulk thermodynamic data as the inputs and give the corresponding GB segregation 
isotherms. 

Figure 5 presents the (relative) chemical potentials and the segregation isotherms for the bulk 
and the exemplar GB with p GB = 0.75 in the Pt-Au system at three different temperatures. The 
chemical potential of the GB is a function of GB density and significantly differs from the 
corresponding bulk chemical potential. The Maxwell construction (defining the two-phase regions) and 
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the spinodal compositions for both GB and bulk materials are built. The graphs show that the GB 
spinodal decomposition can occur for lower values of chemical potential below the bulk spinodal (e.g. 
in Fig. 5 (a)). As a result, the two-phase GB can be in equilibrium with a single-phase bulk. As a matter 
of fact, it has been reported that a segregation-assisted GB phase separation can occur in Pt-Au system 
which results in the formation of two-phase GBs coexisting with a single-phase bulk [68, 69]. 




Figure 5: (a)-(c) The chemical potentials of the Au atoms in bulk and a GB (with p GB = 0.75) Pt-Au are plotted 
for three different temperatures, (d)-(f) Using the parallel tangent condition the segregation isotherms are 
obtained for the corresponding temperatures. Different regimes of segregation appear depending on the GB 
density and temperature. See discussions in Sec. 4.2. 


Furthermore, the results in Figure 5 predicts that the two-phase bulk will be in equilibrium with 
a single-phase GB. In fact, at a given temperature, for any bulk composition within the two-phase region 
of the phase diagram, the bulk is expected to be in equilibrium with a single phase/composition GB. 
These are because of a vertical shift in the GB chemical potential (with respect to the bulk chemical 
potential) that can be seen in Fig. 5 (a). Combining these results with the equilibrium bulk and GB phase 
diagrams (Sec. 3.1), we are able now to study the complete GB thermodynamic phase diagram. This is 
further discussed in Sec. 4. 
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4 . Discussion 

4.1 Coexistence of GB and bulk phases 

In this model, the GB region is described by a continuous atomic density and its spatial gradients. 
The continuity of the atomic density field across GB allows simplifying the Gibbs free energy 
description at the GB plane where the gradient term Vp vanishes (Figure 1 (a) and (b)). The results of 
the current study are obtained for this single material point within the GB region. For the sake of 
simplicity, here we have limited our derivations to the first- and second-order approximations with 
respect to p and Vp, respectively. The current density-based model provides a simple method for 
approximating the Gibbs free energy of a given GB based on available bulk thermodynamic data. The 
model can be applied in full-field simulations of GBs as well, to evaluate the temporal evolution of the 
atomic density and concentration fields. In a parallel study, we have demonstrated the application of 
the density-based model for studying Mn segregation effects in FeMn system [56]. 

In the current study, we focus on the equilibrium GB properties. Using the atomic density-based 
formulation we extend the concept of the phase diagram for GBs, as presented in Fig. 6 (b). We started 
with the equilibrium phase diagram for the bulk and a hypothetical standalone GB, with p GB — 0.75, 
in the Pt-Au solid solution (Figure 4). The results show that like the bulk materials, GBs can also 
undergo a spinodal phase separation that is found to be confined to smaller temperature and composition 
ranges in the Pt-Au system. As the GB density deviates larger from the bulk atomic density, the GB 
two-phase region in the phase diagram becomes smaller. 

The coexistence of the bulk and GB phases can be studied using the equal chemical potential 
condition [20, 67]. Here the bulk phase, which is the dominant part of the microstructure, dictates the 
ultimate thermodynamic states and behavior of the GBs. According to Eq. (22) (or Eq. (23) for the Pt- 
Au system), the GB segregation level depends on the temperature, bulk composition, as well as GB 
type represented by its atomic density p GB . For a regular solution discussed here, the driving forces for 
the segregation can be divided into two parts: The first contribution is the ideal segregation energy 

A E b {^p GB2 — l) + (A K b + pAV B — TAS B ){p GB — 1) that is analogous to the segregation driving 
force in Langmuir-McLean isotherm [67]. If we neglect the energetic contributions due to the mixing 
(fl = 0), we get a density-based version of Langmuir-McLean relation from Eq. (22): 


x gb X B ( AE B (p GB2 - l) + (AK b + pAV B -TAS B )(p GB - 1) 

1-X GB = T^|' eXP l RT 


(24) 


that describes solute segregation for a given GB in an ideal solution. The second contribution in Eq. 
(22) is fl (p Gfi2 — lj + 2[l — p GB2 Xg B | due to the mixing which contains a term with the GB 
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concentration X % B . The two contributions discussed above can have cooperative or competitive effects 
on the segregation level, depending on the magnitude and sign of fl and A E B . 



B spinodally- 

X^ u decomposed GB X^ u 

Figure 6: (a) For Pt-Au system at 700 K, a jump (GB spinodal) occurs before the bulk spinodal decomposition 
takes place. This is a zoomed in plot in Fig. 5 (d) (left comer), see also Fig. 5 (a). (b) Based on the 

current density-based model the coexistence of the bulk and GB phases is depicted in the phase diagram for the 
Pt-Au system. A two-phase bulk will be in equilibrium with a single-phase GB while a two-phase (spinodally- 
decomposed) GB is shown to be in equilibrium with a single-phase bulk. 


Figure 5 presents the chemical potentials and segregation isotherms for the exemplified Pt-Au 
system. In this system, three regimes of bulk-GB coexistence could be identified as a function of the 
bulk composition. First, segregation to the GB occurs as the Au content increases in Pt-rich alloys 
(regime I). For a two-phase bulk, the GB composition becomes fixed (regime II) that applies to the 
entire range of bulk composition in the two-phase bulk region. Any further increase of the Au content 
results in solute depletion in the GB (regime III) with respect to the bulk. In Figs. 5 (a) and (d) different 
regimes of segregation are marked. 

For a dilute Pt-Au alloy (regime I), the GB, which is enriched in Au due to the segregation, can 
undergo a spinodal decomposition before a bulk decomposition becomes possible. At 700 K, for 
instance, a GB spinodal transition is indeed revealed in the Pt-Au system. At this temperature, the GB 
decomposes into low-concentration and high-concertation domains that are in equilibrium with a single¬ 
phase bulk. Figure 6 (a) shows a sharp jump in the GB composition as a function of bulk composition 
(this is a zoomed-in plot of the right corner in Fig. 5 (d)). The jump in the GB composition indicates 
that for a specific bulk composition (in the single-phase region of the bulk phase diagram) a spinodally- 
decomposed two-phase GB will form. Recently, a segregation-assisted GB spinodal decomposition has 
been indeed reported in nanocrystalline Pt-Au alloy [68]. Atomistic simulations of the Pt-Au system 
confirmed the GB spinodal phase transition for bulk compositions well below the bulk spinodal range 
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[69] which is consistent with the current predictions. GB spinodal phase separation is also evidenced in 
other material systems [70] with technological significance [71, 72]. 

For the two-phase Pt-Au bulk materials (in regime II), the current model predicts that a single¬ 
phase GB comes to coexist. This is because of the vertical shift in the GB chemical potential with 
respect to the corresponding bulk chemical potential curve, as depicted in Fig. 5 (a). In practice, a range 
of atomic densities and compositions exists in the GB region (between the GB plane and bulk interior) 
that can enable phase separation in the GB region. Combining the equilibrium GB phase diagram (Fig. 
(4)) and the GB segregation isotherms (Fig. (5)), we can obtain the complete GB phase diagram as 
presented in Fig. 6 (b). 



Figure 7: The GB chemical potential can be ‘Below’ or ‘Above’ the bulk’s curve. Xg S1 and X§ S2 show the bulk 
equilibrium compositions in the two-phase region. The two-phase GB will be in equilibrium with a single-phase 
bulk material as marked by the arrows. The dots on the bulk chemical potential curve indicates the single-phase 

GB is equilibrium with the two-phase bulk. 


The discontinuous jump in the GB segregation isotherm (Fig. 6 (a)) is a result of the difference 
in the free energy functional and chemical potentials of the bulk and GB phases. Depending on the sign 
of A E b and the shape of £l(X B ) function (in composition space), the GB chemical potential and its 
corresponding spinodal decomposition can lie ‘Above’ or ‘Below’ the bulk spinodal (Figure 7). In either 
case, one can activate a segregation-assisted phase separation in the GB region, that can act as a 
precursor for subsequent nucleation, before reaching the bulk two-phase region. Indeed, we have shown 
that a segregation-assisted transient spinodal phase separation can occur at GBs for a larger range of 
bulk compositions [56]. 
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In a polycrystalline material with a large number of GBs of various types (with different 
crystallographic properties), the GB densities are expected to differ from one GB to another. A 
quantitative application of the current model to technical materials, therefore, requires determination of 
the atomic densities for different types of GBs. For this purpose, atomistic simulations can provide the 
necessary means to evaluate GB density profiles. Though atomistic simulation of GBs is beyond the 
scope of this study, in the following we briefly discuss the relationships between the p GB , GB excess 
free volume and GB misorientation angle which can provide some insights to the density-based model. 


4.2 The relation between the GB density and GB nature 

In order to make use of the current model, it is helpful to explore the significance of the average 
GB density and its relation with the GB nature. For a flat GB in a pure substance, the GB width 2ij and 
density p GB relate to the potential energy E A and the gradient coefficient k a . While E A is a property of 
the bulk material, k a maintains the length scale over which the GB is defined. Using atomistic 
simulations, a 3D density map can be extracted from the equilibrium atomic configuration over a 
physical coarse-graining radius that gives 77 and p GB [56]. 

GB energy y is another physical parameter that defines a GB. It is observed that GBs with higher 
energies attract more solute content during the segregation [73]. For a flat GB, Eq. (13) relates the GB 
energy y to its density p GB . Using Eqs. (13) and (18), the GB Gibbs free energy of a regular solution 
can be written as a function of its initial GB energy. Accordingly, we can rewrite Eq. (22) as 


l-X" 


i-xS 


( iAEB + n] (l 1 - i|] 2 - l) + + pAV B - tas b ) ([1 - ig 

■exp - - - 


-1 + 2 n 


'- 11 " t 1 


(25) 


that gives a useful approximation for describing GB’s tendency for segregation as a function of its initial 
energy. 


One way to describe the GB density is to find its relationship with the GB excess free volume 
that is defined as [74] 


AV = 


WN b / 

“wTV 1 


N gb \ 


(26) 
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Here W is the atom diameter, N GB is the number of atoms in the GB, N B is the number of atoms in the 
bulk within the same volume and N 0 is the number of atoms per unit area of the GB. In principle, the 
‘excess’ volume represents the ‘shortage’ of atoms in the GB plane. On the other hand, comparing a 
GB against the bulk of the same volume and using Eq. (12), the relative number of atoms in a GB will 

_ (i+p GB ) 


be n GB = 


(where n B = 1) and the relative difference (with respect to the bulk) will be 


An = 1 — p 


GB 


(27) 


where n — 2 p dr is the (relative) total number of atoms in the GB. Comparing Eqs. (26) and (27) 
one can find the one-to-one relationship between An and AV, where we can identify AV oc ( 1 — p GB ). 
Extensive works have been devoted to calculate and measure GB excess free volume (see for instance 
[75-77]). Aaron and Bolling [74] have shown that the energy of different types of GBs can be associated 
with their excess free volume which can be extended as well for the GB density parameter p GB . 

The lower limit of the GB density value (highest excess free volume) shall be obtained for high 
angle GBs (HAGBs). Arron and Bolling remarked that the excess free volume in a HAGB is comparable 
to that of a liquid phase [74]. On the other hand, the maximum GB density value is the corresponding 
bulk density p B - 1. Special GBs such as coherent twin boundaries (TBs) with low coincidence values 
show small GB energy and excess free volumes [74]. This is well consistent with the current model in 
which as AV goes to zero, p GB approaches p B = 1 and y approaches zero (Eq. (7)). 

For a low-angle GB (LAGB), the situation becomes a bit more complex because of the 
dislocations. Since dislocations are localized defects with a lower atomic density in their core, one 
expects the local GB density p GB within the LAGB plane to fluctuate accordingly. In the absence of 
elastic interactions, we can approximate the average GB density for a LAGB by a simple volume 
averaging of regularly spaced dislocations. For a tilt GB, a simple geometrical construction of a LAGB 

gives the average GB density p GB t ~ 1 — ^2-^ with the GB misorientation 8, and the edge dislocation 
core density p D ~ 0.75 (see Appendix B). For a larger misorientation angle, the GB density will be 
smaller, deviating larger from bulk properties. We also obtain a misorientation-dependent An ~ 
and thus a corresponding GB excess free volume AV oc sin 8 that is indeed consistent with the previous 
studies [74, 78]. GB segregation as a function of GB misorientation has been studied in the Pt-Au 
system. Seki el al. [79, 80] and Seidman [81] studied Au segregation to twist boundaries and have 
shown that although the GB segregation is not homogeneous, its average level increases as the 
misorientation angle increases (with the exception of special GBs). These are consistent with the results 
from our current density-based model where a higher misorientation angle corresponds to a lower GB 
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density. A quantitative study of the GB density values can be conducted using atomistic simulations 
that are left for a future study. 


5. Summary and outlook 

A density-based model for assessing GB thermodynamics is proposed. The current model uses 
the available bulk thermodynamic data as input and gives a rather general formulation for the Gibbs 
free energy of GBs. Using this formulation, the phase diagram of GBs in contact with their 
corresponding bulk phase(s) was obtained. We further discussed the relationships between GB density, 
GB excess free volume and, GB misorientation angle. The model was qualitatively demonstrated in the 
Pt-Au system. The results show that in a Pt-rich alloy Au atoms tend to segregate to the GB. The average 
level of segregation increases with the misorientation angle (with the exception of special GBs) and a 
GB spinodal decomposition can occur for bulk compositions below the bulk spinodal point. A 
comparison to the previous studies on the Pt-Au system confirms these conclusions. 

The current density-based model offers a simple practical approach to develop GB 
thermodynamic databases using the existing bulk thermodynamic data. Despite its simplicity, an 
application of this model can enable knowledge-based studying and design of polycrystalline materials 
by segregation engineering of GBs, as demonstrated in a recent study [56]. A systematic atomistic study 
of the atomic density profiles for different types of GBs can further enrich the current model. Also, the 
current free energy description can be extended by considering elastic energy contribution and a 
generalization to multicomponent alloy systems. 


Appendix A: Equilibrium GB Density Profile 

The equilibrium GB density profile across a GB in a pure substance can be obtained by 
minimizing the GB free energy functional with respect to the relative atomic density. For a ID setup, 


5 $ 

8p 


dp 



= EaP + Cj 


d 2 p 

K w 


= 0 


A.l 


where C/ = + pVj{ — TS% (see Eq. (9)). The general solution for this second-order linear ordinary 

differential equation reads 
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C B 

p (*) = --A + qcos 




where C ± and C 2 are constants and 


dp 

dx 



A.2 


A.3 


Applying the boundary conditions p(x — 0) = p GB . p(x — ±rj) = 1 and the continuity condition 
zr(x = 0) = 0 we can relate the three parameters (E B , C B and K f) ) as 


1 

,-p GB 



2 

A.4 

C 2 = 

0 

A.5 


k a 


u B 

A.6 

1 + P GB 

A.7 


2 



The current derivations allow the continuity of the atomic density profile with Vp(x — 0) = 0. 
This has been also confirmed using atomistic simulations [56]. Common phase-field models for GBs 
with = /[%(1 — (p) 2 + a 2 (V0) 2 ] dV (or other classical phase-field model with O = J[a 1 0 n (l — 
<p) n + a 2 (Vcp) 2 ] dV (n > 1)) which have two free parameters, cannot satisfy the boundary conditions 

<p(x — 0) = Const, <p(x — ±r]) - Const, and the continuity condition = 0) = 0 at the same 
time. Please also see Figure 1 and the discussions in Sec. 2.1. 


Appendix B: A Simple Geometrical Calculation of LAGB Density 

The objective of this appendix is to provide a simple analysis of the LAGB density and its 
dependence on the misorientation angle. We consider the simplest case where the elastic energy of 
dislocations is neglected. A LAGB can be mapped with a set of regularly spaced dislocations separated 
by a distance [63] 
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D - 


sin 9 


B.l 


in which b is the Burgers vector of the dislocations and 8 is the GB misorientation angle. The ‘average’ 
relative GB density can be obtained by averaging the density of a dislocation in the unit area of the 
LAGB. For a tilt LAGB one obtains 


Pmt = 


Db —bb bb r 

-+ — p L 

Db Db H 


B.2 


where p u is the relative atomic density at the core of an edge dislocation. Inserting Eq. B.l in Eq. B.2 


Ptiit = 1 - (1 - P D ) sin 9. 


B.3 


Here it is assumed that the volume between the two dislocations has the same density as bulk ( p B = 1). 


For an edge dislocation depicted in Fig. 8 (a), one can approximate p L 


-(l + —) 

V 4 a-b) 


where a is the 


interatomic distance normal to the edge dislocation line. Neglecting the second term inside the 
parenthesis then we have p D ~ 0.75 that gives the atomic density of a tilt GB as 


Pmt 


1 -■ 


sin 9 


B.4 


In this relation, as 9 approaches zero, p GB approaches p B — 1 and the GB energy (according to 
Eq. (7)) decreases towards zero. Inserting Eq. B.4 into Eq. (13), we obtain 


Ytilt a 0 


(sin oY‘ 
16 


B.5 


that is the tilt GB energy as a function of the misorientation angle. Equations (13), B.4 and B.5 give the 
relationship between the GB density, energy and misorientation angle based on the current simple 

analysis. From Eqs. (27) and B.4, An tat ~ and thus the corresponding average GB excess free 
volume will be proportional to the misorientation angle as AV oc sin 9. Figure 8 (b) and (c) show the 
variation of GB density and energy as a function of the misorientation angle according to this analysis. 

It should be emphasized that Eqs. B.4 and B.5 are based on a very simplified averaging method 
that does not account for the elastic energy of the dislocations which can introduce an error especially 
for obtaining GB energy. Nevertheless, it is useful to show that the current model converges to the 
expected GB energies in the limiting cases, consistent with the previous derivations of GB energy based 
on the excess free volume concept [74]. In order to connect the GB segregation isotherm with its initial 
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misorientation angle, we need to replace the atomic density parameter using B.4. From Eqs. B.4 and 
(19) 


Gg = E B 1 


(l-^) + (K B +pV B -TS B ){l 


sin0\ 


+ 1 


X$ B ^1 (1 - + (Kg + pvg -TSg)(l- 

(>-^i 


£IX BB X BB -TAS^ lx 


Thus, using Eqs. B.6 and (22) the GB segregation isotherm follows 


B.6 


i - x'f 

xl 


l-x; 
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■ exp 
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Figure 8: (a) The core of an edge dislocation has a lower density than that of the bulk that is captured by a 
simple geometrical analysis. GB density (b) and energy (c) as a function of the misorientation angle are plotted 
for a tilt LAGB (Eqs. B.4 and B.5) with the exception of special GBs. 
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